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1.0  INTRODUCTION 


1.1  Purpose 

This  report  presents  the  theory  of  the  dynamic  response  of  a  linearly  elastic, 
circular  arch,  including  the  effects  of  both  shear  deformation  and  rotary 
inertia.  The  arch  central  angle  is  arbitrary,  and  the  end  restraints  are 
linearly  elastic.  The  theory  is  extended  to  include  linear  viscoelastic  beha¬ 
vior,  in  which  both  arch  and  end  restraints  are  governed  by  the  same  single 
linear  rate  sensitivity  parameter. 


1.2  Motivation 

The  work  described  herein  is  an  analytical  extension  of  results  presented  in 
Reference  1.  The  approach  was  motivated  by  Dr.  Timothy  Ross'  success  in 
explaining  direct  shear  failure  of  dynamically  loaded  reinforced  concrete  box 
roof  slabs  using  a  Timoshenko  beam  model  and  an  appropriate  reinforced 
concrete  direct  shear  failure  criterion  (Ref.  2). 


1.3  Application 

When  the  theoretical  results  presented  herein  have  been  evaluated  numerically 
for  prescribed  arch  central  angle,  end  restraint,  and  loading,  they  will  apply 
to  the  KACHINA  arches  tested  at  the  Air  Force  Weapons  Laboratory  (Refs.  3-5). 
Admittedly  reinforced  concrete  is  not  always  linearly  elastic,  so  application 
of  the  theory  to  the  KACHINA  arches  will  require  judgment.  Nevertheless,  the 
convenience  of  a  closed  form  solution  afforded  by  elastic  theory  is  too 
attractive?  not  to  exploit  for  the  perspective  it  gives  on  overall  dynamic 
behavior  and  the  influence  of  structural  and  loading  parameters. 

1.4  Point  of  Departure 


The  point  of  departure  for  the  dynamic  analysis  is  the  set  of  equations  for 
the  static  behavior  of  a  circular  planar  member  presented  by  Connor  in  Chap¬ 
ter  14  of  Reference  6.  Other  related  work  is  presented  in  References  7-11. 


Details  of  the  dynamic  .analysis  are  contained  in  Appendices  A-L.  The 
discussion  in  Section  2  below  summarizes  each  appendix  in  order.  Because  of 
the  detailed  treatment  and  the  many  equations,  a  few  letter  symbols  have  been 
used  differently  in  separate  appendices.  However,  symbols  are  defined  where 
first  used  in  each  appendix,  so  there  should  be  no  confusion.  See  Appendix  M. 

2.0  DISCUSSION 

Figure  1  shows  a  portion  of  a  circular  arch,  and  defines  the  terms  which 
control  its  dynamic  response.  The  arch  centerline  radius  is  R,  and  the 
central  angle  defining  a  particular  plane  normal  cross  section,  PQ,  in  the 
undeformed  arch  is  8.  The  width  and  depth  of  the  constant  rectangular  arch 
cross  section  are  b  and  d.  Deformation  of  the  arch  causes  the  centroid  of  a 
plane  cross  section  PQ  to  displace  tangentially  by  an  amount  u^  and  radially 
inward  by  an  amount  U2,  and  the  cross  section  to  rotate  through  an  angle  <p; 
but  the  section  is  assumed  to  remain  plane,  so  that  P*Q*  is  a  straight  line. 
The  distributed  tangential,  radial,  and  moment  loads  acting  at  point  0  in  the 
undeformed  arch  are  b^,  b2,  and  m.  Notice  that  the  angle  is  an  independent 
displacement  parameter,  so  that  P*Q*  is  generally  not  normal  to  the  deformed 
centroidal  axis.  The  angular  amount  by  which  section  P*Q*  deviates  from  being 
normal  to  the  deformed  centroidal  axis  is  the  arch  shear  deformation  at 
point  0. 

Appendix  A  begins  by  formulating  the  equations  of  motion  for  a  differential 
arch  element,  using  the  internal  stress  resultants  (Equations  A14-A16).  Next, 
the  form  of  the  displacement  field  is  defined  by  assuming  that  plane  cross 
sections  normal  to  the  undeformed  arch  centroidal  axis  remain  plane  during 
deformation.  This  is  equivalent  to  expanding  the  displacements  in  a  Taylor 
series  in  two  variables  about  the  centroid  (arc  length  and  inward  radial  posi¬ 
tion)  and  retaining  only  the  constant  and  linear  terms.  Having  assumed  the 
form  of  the  displacement  field  (Equation  A23),  the  extensional  and  shear 
strains  are  obtained  as  functions  of  the  displacement  field  parameters  and 
their  derivatives  with  respect  to  arc  length  (Equations  A31-A33).  Three 


strain  parameters  are  defined  at  the  centroid  (Equations  A34-A36),  and  then 
the  strain  components  anywhere  in  the  arch  are  expressed  as  functions  of  the 
strain  parameters  (Equations  A37-A39).  Normal  and  shear  stresses  are  obtained 
as  functions  of  the  strain  parameters,  using  elastic  stress-strain  equations 
(Equations  A40  and  A41),  and  then  integrated  over  the  cross  section  to  obtain 
the  internal  stress  resultants  as  functions  of  the  strain  parameters 
(Equations  A42-A44).  The  stress  resultant  equations  of  motion,  stress 
resultant-strain  component  equations,  and  strain* component-displacement 
equations  are  summarized  in  matrix  form  in  Equations  A61,  A63,  and  A64.  Use 
of  matrix  notation  greatly  simplifies  the  entire  analysis.  It  is  compact, 
helps  avoid  algebraic  errors,  and  even  suggests  analytical  approaches  not 
obvious  when  equations  are  written  out  in  detail. 

The  equations  of  motion  in  terms  of  displacement  components  and  their  deriva¬ 
tives  are  obtained  in  matrix  form  in  Appendix  B,  by  substituting  Equation  A64 
into  Equation  A63  and  that  result  into  Equation  A61.  The  final  result  is 
Equation  B17. 

Appendix  C  considers  the  homogeneous  (free  vibration)  form  of  the  equations  of 
motion  (C2),  and  uses  the  classic  separation  of  variables  approach  to  define  a 
vibration  mode  (Equation  C3).  Free  vibration  in  a  single  mode  is  shown  to  be 
harmonic  (Equation  C15),  and  the  mode  shapes  are  shown  to  obey  a  coupled  set 
of  three,  linear,  ordinary,  second  order  differential  equations  (CIS) .  The 
three,  coupled,  mode  shape  differential  equations  are  uncoupled  by  pre¬ 
multiplying  them  by  the  adjoint  of  the  coefficient  operator  matrix  (Ref.  12), 
but  at  the  cost  of  having  the  resulting  three  separate  identical  linear  ordi¬ 
nary  differential  equations  be  of  order  six  instead  of  two  (Equations 
C21-C23) . 

Equation  C20  defines  the  sixth  order,  linear,  ordinary  differential  operator, 
A,  which  appears  in  each  of  Equations  C21-C23.  These  are  the  spatial  dif¬ 
ferential  equations  which  the  mode  shapes  must  satisfy.  The  operator.  A,  is 
defined  as  a  determinant,  which  when  expanded  yields  a  sixth  order  polynomial 
in  the  frequency  parameter,  a,  (Equation  D17)  or  a  sixth  order  polynomial  in 
the  spatial  differential  operator,  0,  (Equation  D18).  '  Equation  D18  is  the 
more  fundamental  form,  and  shows  the  frequency  dependence  of  its  coefficients. 


The  last  half  of  Appendix  D  discusses  the  way  in  which  the  coefficients  in 
Equation  018  vary  with  frequency. 

The  operator,  A,  defined  by  Equation  C20  is  a  sixth  order,  linear,  ordinary 
differential  operator,  but  the  odd  order  derivatives  are  missing.  Therefore, 
assuming  the  solution  of  the  homogeneous  equation  AX  =  0  (Equation  C21)  to  be 
of  the  form  X  =  yields  a  cubic  character istic  equation  in  X^  (E6). 

Whether  the  parameter  X  is  real,  imaginary,  or  complex  depends  on  the  co¬ 
efficients  of  the  cubic  equation.  Appendix  E  discusses  the  algebraic  solution 
of  a  cubic  equation,  as  the  basis  for  a  detailed  examination  of  the  ncture  of 
the  characteristic  roots  and  their  associated  mode  shape  functions,  which 
appears  in  Appendix  F. 

Appendix  F  examines  the  fifteen  possible  characteristic  root  combinations  and 
their  associated  mode  functions.  The  mode  functions  are  a  closed  set  of  six 
linearly  independent  functions,  in  the  sense  that  differentiation  of  any  one 
yields  a  linear  combination  of  the  other  five. 

Appendix  G  establishes  the  derivative  properties  of  symmetric  and  antisym¬ 
metric  functions  needed  to  define  arch  centerline  boundary  conditions  for 
symmetric  and  antisymmetric  modes. 

Appendix  H  explains  in  detail  how  the  natural  frequencies  and  associated  mode 
shapes  are  calculated  for  given  arch  boundary  conditions.  The  feature  which 
distinguishes  this  modal  analysis  from  many  others  is  the  fact  that  the  pre¬ 
cise  form  of  the  frequency  equation  is  frequency  dependent. 

Appendix  I  presents  the  arch  modal  orthogonality  relations.  The  approach  is 
that  for  a  Sturm-Liouvi lie  problem,  supplemented  by  Rayleigh's  ingenious 
application  of  L'Hospital's  Rule  to  find  the  integral  of  the  inner  product  of 
a  mode  with  itself  over  the  arch  length  (Refs.  13  and  14). 

Once  the  free  vibration  mode  orthogonality  relations  have  been  established, 
the  formal  solution  of  a  transient  forced  vibration  problem  becomes  straight¬ 
forward.  Both  the  arch  displacements  and  distributed  loads  are  assumed  to  be 
expressible  as  modal  series  expansions,  and  the  equations  of  motion  are  then 
scanned  with  a  particular  mode  shape.  The  result  is  a  single  degree  of 
freedom  differential  equation  for  the  associated  modal  amplitude  (J6),  the 
solution  of  which  is  a  Guhamel  convolution  integral  (Equation  J8). 


Appendix  K  extends  the  above  elastic  arch  analysis  to  the  case  of  viscoelastic 
behavior.  When  both  arch  and  boundary  restraints  are  governed  by  the  same 
single  rate  sensitivity  parameter,  the  elastic  modal  analysis  still  applies. 


Appendix  L  addresses  the  most  difficult  computational  phase  of  arch  modal 
analysis,  finding  the  roots  of  the  frequency  equation  (H7).  The  frequency 
equation  is  a  complicated  transcendental  equation,  usually  involving  both 
trigonometric  and  hyperbolic  terms.  The  modal  frequency  parameter,  o,  not 
only  appears  implicitly  in  several  places,  but  also  controls  the  precise  form 
of  the  frequency  equation.  Considerable  care  is  needed  in  calculating  the 
higher  modal  frequencies,  because  some  of  them  are  apt  to  be  closely  spaced. 


3.0  CONCLUSIONS 

The  distinguishing  feature  of  the  closed  form,  dynamic  arch  analysis  presented 
in  this  report  is  its  extensive  use  of  matrix  notation.  The  analysis  includes 
the  effects  of  both  shear  deformation  and  rotary  inertia  on  the  transient 
response  of  both  an  elastic  and  a  viscoelastic  circular  arch,  having  an 
arbitrary  central  angle  and  corresponding  elastic  or  viscoelastic  end 
restraints.  The  resulting  equations  are  presented  at  a  level  of  detail  suf¬ 
ficient  for  direct  computer  programming,  and  a  computer  program  (a  modifica¬ 
tion  of  ZEROIN)  is  presented  for  finding  the  modal  frequencies. 
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APPENDIX  A 


GOVERNING  EQUATIONS  FOR  THE  DYNAMIC  RESPONSE  OF  A 
CIRCULAR  ARCH,  INCLUDING  BOTH  SHEAR  AND  ROTARY  INERTIA 


Figure  A1  shows  a  differential  element  of  a  circular  arch.  For  the  sign 
convention  shown,  with  right-handed,  orthogonal  unit  vectors  t  •,  and  t^. 


t3  =  t  x  tj  =  constant 


[Al] 


1 

R 


[A2] 


If  we  set 


F  =  Fjtj 

+  F2*2 

(internal  force  resultant) 

[A4] 

(internal  moment  resultant) 

[A5] 

b  =  b1T1 

+  b2*2  . 

(external  distributed  load) 

[A6] 

m  =  mt3 

(external  distributed  moment) 

[at] 

“  ■  Vi 

+  V2 

(centroidal  displacement) 

[A8] 

5  =  *t. 

(cross-section  rotation) 

[A9] 

then  the  translational  equation  of  motion  for  the  element  is 


F  + 


r  3F  ^ 

F  +  —  ds 
3s 


+  bds  =  pAdsu 


[A10] 


where  p  is  mass  density,  and  A  is  cross-sectional  area. 


.WWV  V  11'  W'V.l'WT.'  W  W^.V.7.7JVIV.’I  V.V.V.V.V.’V  7  V  ■  V  V  ■ 


or 


3F  - 

-  +  b  =  pAu 


and  the  rotational  equation  of  motion  for  the  element  is 


-  M  + 


-  3M 
M  +  —  ds 
3s 


+  ds  x  F  +  mds  =  plds4» 


where  I  is  the  moment  of  inertia.  Then 


3M 


—  +  t,  x  F  +  m  =  pI4> 

oS  1 


Substitution  of  Equations  A2  through  A9  into  Equations  All  and  A13 
yields 


3F1  F2 

3s  R  +  bl  pAui 


3F2  F1 

aT  +  r  +  b2  s  pAu2 


3M 


3s  *  F2  *  m  *  P I* 


[All] 


[  A 1 2  ] 


[  A 1 3  ] 


[A14] 


[A15] 


[a  16] 


.The  next  step  in  deriving  the  governing  equations  is  to  obtain  expressions 
for  the  longitudinal  extensional  and  shear  strains  at  any  point  in  an  ori¬ 
ginally  normal  cross  section.  Consider  a  point  Q  in  the  undeformed  arch 


shown  in  Figure  A2,  the  position  vector  for  which  is  denoted  r^  in  Figure 


A3.  From  Figures  A2  and  A3  we  have 


r,  =  r  >  y2t2  ♦  y3t3 


[ait] 


so  that  in  Figure  A3, 


PP.  =  r1  ds  = 


t.  ds 


fulfil 
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The  centroidal  displacement  vector  of  point  P  in  the  deformed  arch  is 


“  ■  Vi  *  V2 


and  the  assumption  that  plane  cross  sections  normal  to  the  centroidal  axis 
before  deformation  remain  plane,  and  normal  to  the  1-2  plane  during  de¬ 
formation  yields  an  expression  for  the  displacement  vector  of  point  Q  in 
the  deformed  arch 


u2  =  u  +  x  y2t2  =  (u1  -  <i^2)t1  +  u2t2 


*  *  *  * 

If  P  ,  Q  ,  Q  ,  and  Q  denote  the  displaced  points  P,  Q,  Q  ,  and  Q  in  the 
deformed  arch,  then  the  position  vectors  of  points  P  and  Q  are 


-*  -  - 
r  =  r  +  u 


r2  *  r2  +  U2 


so  that  in  the  deformed  arch.  Equations  A25  and  A19  yield 


[A25] 


[•Ml] 


'A  ^  ^ sV  s*  v  W.v 


and  Equations  A25  and  A20  yield 
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Assuming  the  displacement  derivatives  are  small  with  respect  to  unity 
permits  neglect  of  the  nonlinear  terms  in  Equations  A28,  A29,  and  A30, 
which  leaves,  using  Equations  A21  and  A23, 
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Note  that  y  is  positive  when  the  dot  product  Q  Q  •  Q  Q  is  positive, 

^  *  if  if  i  C 

which  means  that  the  angle  Q^Q  Q ^  is  less  than  n/2.  When  this  is  the 
case,  the  direction  of  the  shear  stresses  on  a  deformed  arch  element  are  as 
Shown  in  Figure  A4 . 

If  we  now  set 


1  2 

=  (e  1  =  -  -  — 

1  y2=0  3s  R 

[  A3  4  ] 

U1  3u2 

=  ,7u)y2-o  -  -  *  *  r  *  5 r 

[  A3  5] 

djf 

3s 

[A36] 

then  Equations  A31,  A32,  and  A33  can  be  written  in  the  form 


Assuming  zero  initial  strain,  the  longitudinal  normal  and  transverse  shear 
stresses  are 
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11 


y2k) 
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12 


[A40] 


[A41] 


where  E  is  Young's  modulus,  and  G  is  the  shear  modulus,  so  that  the  stress 
resultants,  F^ ,  F 2,  and  M  in  Figure  A1  are 


fo,udy2dy3 


=  E6 


rr_^2 


UJ 


y,  dy3 


Ek 


F„  = 


ff 


ff_^2 


JJa12dy2dy3  =  G62jj - 7  dy3 

1  -  r 


M  3  -\\ 


c^y^y^y,  =  _E6 
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y,  dy3 
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The  three  integrals  appearing  in  Equations  A42,  A43,  and  A44  can  be 
expressed  in  terms  of  a  single  integral.  Referring  to  Figure  A2,  we  write 


2J  1  -r 


fl  +  xl  R  .  fl  ♦  xl 

,n[r^j =  [a  1nlr^]|< 


.  fl  ±  xl  . 

=  lni“J  .  .  fd! 


=  [g-]d  =  Cjc 


[A51] 


where 


41  -* :} 


2  4  6 

.  XXX 

1  +  r  *  r  +  r  * 


[A52] 


and  therefore 


(Vd/R=0  =  1 


Equation  A47  can  be  written  in  the  form 


■ 

* 

y2dy2 

5  -1 

Y2 

d 

1  •  r 

R 

d2  =  C„d2 


[A54] 


where 


and  therefore. 


(Vd/R=0  =  0 


Equation  A48  can  be  written  in  the  form 


/  -•  /  -•  .  .*  . 


Vi  X  x!  x! 

2x  =  6  +  10  +  14  + 


[A55] 
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and  therefore. 


(C3}d/R=0  =  1 


[A59] 


Values  of  C^,  C£,  and  C3  are  tabulated  below  to  four  significant  figures  as 
a  function  of  d/R. 


d/R 

C 

1 

C 

2 

C 

3 

0 

1.000 

0 

1.000 

.10 

1.001 

.0083 

1.002 

.20 

1.003 

.0168 

1.006 

.30 

1.008 

.0253 

1.014 

.40 

1.014 

.0342 

1.025 

.50 

1.022 

.0433 

1.039 

If  we  define  the  angular  differential  operator,  D,  to  be 


2-  .  0 

30 


then  Equations  A14,  A15,  and  A16  can  be  written  in  the  form 


36 


Equations  A42,  A43,  and  A44  can  be  written  in  the  form 


where  v  is  Poisson's  ratio.  Equations  A34,  A35,  and  A36  can  be  written  in 
the  form 


Equations  A61,  A63,  and  A64  are  the  governing  equations  for  the  dynamic 
response  of  a  circular  arch,  including  both  shear  and  rotary  inertia,  with 
one  exception.  That  exception  is  replacement  of  the  constant  Cj  in  the 
second  of  Equations  A63  by  the  transverse  shear  coefficient,  k',  to  account 
for  the  fact  that  the  arch  cross  section  does  not  remain  plane  (Ref.  15). 


Ref.  15.  Oden,  J.  T. 

2nd  Edition 


and  E.  A.  Ripperger, 
McGraw-Hi 11,  ( 1981 ) 
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The  equations  of  motion  (A61)  then  take  the  form 
J  (01  ♦  G ) { F }  +  {b}  =  pA{u) 

where  X  is  the  identity  matrix, 

the  stress  resultant-strain  equations  (A63)  take  the  form 
{F|  =  EAS  {6 |  . 

and  the  strain-displacement  equations  (A64)  take  the  form 
{6}  ■  J  (DI  -  G T ) { u } 

Substitution  of  Equations  B8  and  B9  into  Equation  B7  yields 

^  (DI  +  G)S(DI  -  G T ) { u |  +  (b|  =  pA  fu) 

R^ 

or 

(DS  +  GS)(DI  -  GT){u|  +  (bf  =  pA  {u| 

R 

or 

^  02S  -  0SGT  +  OGS  -  GSGT  {uj  +  (b|  =  pA  {uj 

R 

or 

02S  +  0(GS  -  SGT)  -  GSGT  {u}  +  {b }  =  pA  {u| 

R^ 


The  matrices  in  Equation  Bll  are  computed  below. 


2(1  ♦  v) 


k' 

2(1  ♦  v) 


-C2  f 


k'R 


2r( 1  +  v) 


k'R 

2r( 1  ♦  v) 


2r(l  +  y) 


2r  ( 1  +  v) 


If  we  set 


GSG  =  K 


then  Equation  Bll  can  be  written  in  the  form 


{u}  -  [jj]  (02S  +  OH  -  K){u|  *  ^  {b} 


[B17 


Equation  B17  is  the  governing  equation  for  the  dynamic  response  of  a 
circular  arch,  including  both  shear  and  rotary  inertia,  written  in  matrix 
form.  The  matrices  S  and  K  are  symmetric,  and  the  matrix  H  is  antisym¬ 
metric  ,  i .e. , 


ST  *  S 


«T  =>  K 


HT  =  -H 


Note  also  that  Equation  A62  yields 


S  =  2/3 
r 


c^a.v.v.s 


APPENDIX  C 


UNCOUPLING  THE  EQUATIONS  OF  MOTION 


For  free  vibration,  we  set 
{b}  =  {0} 


in  Equation  BIT,  and  obtain 


r  I2  ? 

{u}  -  §  (0  5  +  OH  -  K) { u }  =  {0} 


Assuming  a  separable  solution,  we  set 


M  =  T{U| 


where  T  is  a  function  of  time  only,  and 


{U}  =  Y 


where  X,  Y,  and  Z  are  functions  of  0  only. 


Substitution  of  Equation  C3  into  Equation  C2  yields 


Tiui  -  f|!  T(02S  ♦ 


OH  -  K)(U)  =  |0( 


Cel2 

T{U}  -  [pj  TL(U')  «  (0) 


where  L  is  the  linear,  spatial,  matrix  differential  operator, 


L  *  OS  ♦  OH  -  K 


It  is  convenient  to  define  a  diagonal  matrix  containing  the  elements  X,  Y 


V  TTTT^TT^yrTr 


r. 

e. 

« 

k, 

& 

£ 

§ 


! 


,1 


i 


'  % 


w. 

& 


r" 

£ 


fa,  - 


X  0  0 
0  Y  0 
0  0  Z 


(The  notation  f  j  denotes  a  diagonal  matrix.) 
as  well  as  the  identity  vector 


{m}  = 


so  that 


Tm 


-1 


and 


Uj  (U}  =  jm} 


rujH  =  |u| 


f-  _1 

Premultiplication  of  Equation  C6  by  'Uj  yields 


T{m|  - 


2  r-  _1 

T'Uj  L(U}  =  {0} 


and  division  by  T  yields 


fM  ■ 


2r- _l 

’u.  L{U} 


[  C8  ] 


[C9] 


[CIO] 


[Cll] 


[  C 1 2  ] 


[  C 1 3  ] 


Since  the  LHS  of  Equation  C13  is  a  function  of  time  only,  and  the  RHS  is 
a  function  of  0  only,  for  arbitrary  values  of  T  and  0,  both  sides  must  be 
constant,  so  that 


f{m|  = 


2r-  'l  ? 

U  j  Lju}  =  V(m} 


[C14] 


•;.V 

•m\A 
•>  •*» 

i 


>v-: 

:g| 

‘.V/1 


r.  ■r.' 

W 

'C<-' 


.S  V 
V  J 

I 

.-.V 

V 
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*> 

r*1 


y. 


a*  ^  -  -  • 


and  therefore 


T  +  p  T  =  0 


The  function  T  is  therefore  harmonic.  Premultiplication  of  the  last  two 
parts  of  Equation  C14  by  ^Ui  now  yields 


|  klu}  =  -p'{u} 


where 


L{U}  =  -  -^{U}  = 


-c2{u} 


Equation  C16  can  be  written  in  the  form 


(L  ♦  aZI)|U}  =  {0} 


The  linear  differential  operator  coefficient  matrix  (L  +  ocl)  in  Equation 
C18  is  coupled,  but  the  system  of  ordinary  differential  equations  can  be 
uncoupled  by  premultiplying  by  the  adjoint  of  the  coefficient  matrix. 

(L  +  *  <72I)|U}  =  Alju}  =  A jl)}  =  jo}  [C19 

where  A  is  the  determinant  of  the  coefficient  matrix  in  Equation  C18. 


A  =  l  ♦  a  I 


W>  . >  ».  *  ■>  fa  •  O 


L-w  ’.>Vn^  vv v  • 


APPENOIX  D 


EXPANSION  OF  THE  CHARACTERISTIC  DETERMINANT 


In  Equations  C21,  C22,  and  C23,  the  linear,  ordinary,  differential 
operator  A  represents  the  expansion  of  a  determinant. 


A  =  L  +  <r  I 


[C20  bis] 


where 


L  *  0  S  ♦  DH  -  K 


[C7  bis] 


-2K3C, 


-2/3C, 


k' 

2(1  +  y) 


2(1  +  u) 


k' 

2(1  +  v)  *  C1 


k'R 

2r ( 1  +  y) 


•W.01 

.V.V’I 


■■■ 1  +  2/3C 
2(1  +  y)  2 


k  * 

2(1  +  y) 


k'R 

2r ( 1  *  y) 


[B15  bis] 


k'R 

2  r ( 1  ♦  y] 


To  facilitate  the  expansion  of  A,  let 


C1  =  A 


2r2 ( 1  +  y)  j 


-i 

sX'/A 


vy:< 

-.-■.'■J 
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2(1  +  tf) 


C3  =C 


2/3C2  =  e 


5  =  F 
r 


Then  Equations  C7,  Dl,  D2,  and  B15  yield 


L  =  D 


A 

0 

-E 

0 

-(A+B) 

0 

0 

B 

0  +  D 

(A+B) 

0 

-(BF+E) 

-E 

0 

C 

0 

(BF+E) 

0 

B 

0 

-BF 

0 

A 

0 

-BF 

0 

BF2 

ad2-b 

-  (A+B)0 

-ed2+bf 

(A+B)0 

< 

i 

cv 

o 

CO 

-{BF+E)D 

-ED2+BF 

(8F+E)D 

2  2 
CD^-BF^ 

so  that  Equation  C20  takes  the  form 

I  (A02-B)+o2  -(A+B)D 


(A+B)0 


-ED  +BF 


2  2 
(BO  -A)+<7 


(BF+E)D 


-ED  +BF 

- (BF+E ) 0 

2  2  2 
(CD  -BF  )+o 


Equation  D9  is  very  similar  to  the  determinant  which  arises  in  calculating 
principal  stresses,  so  that  the  expansion  of  Equation  D9  is  known  to  take 
the  form 


A  =  o6  ♦  la4  +  I.ct2  +  I. 


roioi 


where 


Ij  =  Tr  (L) 


[Dll] 


2 

2 

AO  -B 

-(A+B)0 

A 

BD  -A  -(BF-t-E)D 

(A+B)0 

bo2-a 

▼ 

(BF+E)D  CD2-BF2 

2  2  2 

CD  -BF  -ED  +BF 

-£02+BF  A02-B 


ad2-b 

(A+B)0 

-ED2+BF 


-(A+B)D 

2 

BD  -A 
(BF+E)D 


-ED  +BF 


-(BF+E)D 

2  2 
CD  -BF 


[012] 


[013] 


Equations  Dll,  D12,  and  D13  yield 

Ij  *  ( A  ♦  B  C2  )  D  -  (A  *  B  +  BF2)  [014] 

12  (AB  +  BC  ♦  AC  -  E2)04 

♦  ( 2AB  -  AC  ♦  4BEF  +  E2  -  BC  -  ABF2)02  ♦  (AB  ♦  ABF2)  [D15] 

13  =  (ABC  -  BE2 ) (O2  «•  1)202  [016] 

Substitution  of  Equations  D14,  D15,  and  D16  into  Equation  DIO  yields 


A  =  a6  +  [(A  +  B  +  C)D2  -  (A  +  B  ♦  BF2)]ct4 

+  [(AB  +  BC  +  AC  -  E2)D4  +  (2AB  -  AC  +  4BEF  +  E2  -  BC  -  ABF2)D2 
+  (AB  +  ABF2)](72 

+  [(ABC  -  BE2)06  +  2 (ABC  -  BE2)04  +  (ABC  -  BE2)D2]  [D17] 


Finally,  A  can  be  expressed  as  a  polynomial  in  the  operator  D. 


A  =  (ABC  -  BE 


2)06  +  [(AB  +  BC  + 


AC  -  E2)a2  +  2 (ABC  -  BE2)]d4 


+  [(A  +  B  +  C)a4  +  (2AB  -  AC  ♦  4BEF  ♦  E2  -  BC  -  ABF2)o2 
+  (ABC  -  BE2 ) ] D2  +  [a6  -  (A  +  B  +  BF2)o4  +  (AB  ♦  ABF2)a2]  [D18j 


All  the  coefficients  in  Equation  018  can  be  expressed  in  terms  of  the 

t  h 

invariants  of  the  matrices  S,  H,  K,  and  S  -  K,  where  I ^  (x)  is  the  j 
invariant  of  the  matrix,  x. 

I.  (S)  »  A  ♦  B  ♦  C  | 

I2  (S)  =  AB  ♦  BC  +  AC  -  E2  | 

I  (S)  =  ABC  -  BE2  | 

o 

Ij  (H)  =  0  j 

12  (H)  =  A2  +  2A8  +  82  +  82F2  +  2BEF  +  E2  | 

13  (H)  «  0  | 


ll  (K)  =  A  +  B  +  BF 


I2  (K)  =  AB  ♦  A8F* 

r3  <*>  "  0 


[019] 

<  f  » 

[D20] 

*.v\ 

[D21] 

'S. 

[D22] 

V* 

$S: 

[D23] 

V/A 

.‘aa 

[024] 

■  v.v. 

[025] 

.'AA 

.  A 

AA 

ir 

[026] 

& 

•  V  */ 

[027] 

A-8 

S  -  K  =  0 


I2  (§  -  K)  =  -A2  +  2AB  -  B2  -  B2F2  +  2BEF 


[D29] 


I,  (S  -  K)  -  I2(S)  -  I?(K)  +  I  (H) 


2A8  -  AC  +  4BEF  ♦  E  2-  BC  -  ABF  2 


[D30] 


-'-A  v'.-'  - •  -■  .  -  .  ■  ■ v  ■  ■ 


Thus,  Equation  018  can  be  written  in  the  form 


A  =  I3  (S)D6  +  I2(S)a2  +  2I3(S)  D4 


+  [I1(S)0  +  [  I2(S  -  K)  -  I2(S)  -  I2(K)  +  I2(H)]aZ  ♦  13(S)Jd': 

+  ct6  -  ^(KJa4  +  I2(K)a2  [D31 


Substitution  of  Equations  D3  -  07  into  Equations  D19,  D20,  021,  025,  D26, 
and  030  yields 


h  '§>  ■  is  *  S|  *  Trrrn 
S  's'  *  irrhy  [s  *  s]  .*  [ss  -  12c2 
S  <®>  ■  JiTT^T  [ss  -  12C2 

S  *  s  *  Jirbr  [‘  *["f. 


s  <5»  -  ittVitt  i  * " 


2  (s  -  K)  -  r2(S)  -  i2(K)  +  r2  (H) 


2TT77)  [2Ci  -  S  ♦  ^c2[*]  -  ci[" 


-  ci 


C3  ~  12C2 


Therefore,  if  we  set 


V?)  2 

a  =  vi)a  +2 


I,(S)0  +  fl0(S  -  K)  -  T  (S)  -  I„(K)  +  U(H)1ct" 


Using  the  commonly  accepted  value 


k.  _  5  .si 

*  "  6  ~  12 


Equation  051  yields 

„  _  12(1  +  jj] 


u 

n 

0.10 

2.64 

0.15 

2.76 

0.20 

2.88 

0.25 

3.00 

0.30 

3.12 

0.35 

_ 

3.24 

Assuming,  based  on  the  values  tabulated  on  the  previous  page,  that 
K0  =  K1  =  K3  3  1 

k2  ■  0 


Equations  053,  054,  and  055  reduce  to 


a  =  (2  +  n)<j  ♦  2 


b  *  (2n  +  l)a4  -  (z2  ♦  n  -  l)c2  +  1 
c  =  [no4  -  (z2  +  n  +  l)o2  +  (1  +  z2)]a2 


and  if  we  assume  that 
n  *  3 


2  2 
z  ♦  n  ±  1  =  z 


Equations  058-060  reduce  to 


a  =  5a  ♦  2 


APPENDIX  E 


GENERAL  SOLUTION  OF  THE  CHARACTERISTIC  DIFFERENTIAL  EQUATION 


The  constant  I3  (S)  in  Equation  041  is  nonzero,  and  can  therefore  be 
ignored.  We  therefore  seek  the  general  solution  of  the  homogeneous, 
linear,  ordinary  differential  equation 


-X  =  0 


where 


L  =  06  +  aO4  +  bO2  +  c 


0  =  h 


The  solution  of  Equation  El  is  assumed  to  be  of  the  form 
xe 


X  =  Ke 

and  substitution  of  Equation  E4  into  Equation  El  yields 


(X6  +  aX4  +  bX2  +  c)KeX0  =  0 


If  a  nontrivial  solution  exists,  it  must  be  that 


,6  '  4  2 

X  +  aX  +bX  +c=0 


Equation  E6  is  a  cubic  equation  in  X2,  and  therefore  has  three  roots 


2  2  2 
X j ,  X 1 1 1  X 1 1 1 


Equation  E6  can  therefore  be  written  in  the  form 


[x2-A*lp-A 


,2  ,2  *6  .4  2 

X  -  *■  +aX  +bX  +  c  =  0 


Comparison  of  Equations  E6  and  E7  shows  that 


a  =  - 


2  2  2 
XI  +  XII  +  XIII 


b  = 


,2.2  ,2  .2  .2  2 

III  +  XII  XIII  +  XIIIXI 


.2  .2  .2 

C  "  "  I  II  III 
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then  Equations  E16  and  E17  will  be  satisfied  by  setting 


y  .  =  A  .  +  B  . 
J  J  J 


If,  on  the  other  hand, 


(j  =  1.  2, 


3) 


[ 


[ 


then  the  quantity  R  defined  by  Equation  E20  is  negative,  so  that  /r  is 
3  3 

imaginary  and  A  and  B  are  complex.  Then  if  we  set 


3  3 

then  we  can  show  the  complex  quantities  A  and  B  in  an  Argand  diagram 
(Figure  E).  Note  that  Equation  E33  guarantees  that  P  is  positive. 
Figure  E  yields 


where 


cos  3u  = 


9 

2 


3/2 


If  we  set 


Wj  =  u 

2rr 

(ji  =  M  -  - 

2  3 


27T 

u„  =  u  +  — ■ 
3  3 


and  then  set 


A  . 
J 


bks 


1/2  iu. 
e  J 


1/2  iu. 


(0  <  3u  <  n) 


[E38] 


[E39] 

[E40] 

[E41] 


(j  =  1.  2,  3) 


(j  =  1.  2,  3) 


then  Equations  E16  and  E17  will  be  satisfied  by  setting 


y .  =  A .  +  B .  =  2 J  -  cos  u  . 
J  J  J  '  3  ■  J 


(j  =  1.  2,  3) 


t  E  4-2  ] 


[E43] 


[  E  4  4  ] 


I* 


-  " 


i 


mi 


.•-»j 

■lyA: 

» v  v, 

•W, 

V  Ar 

•>>o 

brd 


.■ '/ V 


Finally,  in  keeping  with  Equation  Ell,  we  write 

,2  a 

XI  "  Y1  '  3 

[E45a] 

.  %  \  i 

-J 

.2  a 

XII  ”  y2  "  3 

[E45b] 

N.", 

-H  • •  o 

.2  a 

XI I I  ~  y3  "  3 

[E45c ] 

\*  J 
»^VJ 

2  2  2 

We  now  need  to  determine  whether  the  roots  ,  *11'  anci  *IH  are  rea^  or 
complex,  and  if  real,  whether  they  are  positive  or  negative. 


:-V\ 

v\ 

'.y/j 


*.  •.  v 

y;-\ 


p  .*  > 

/>\v 


.  v  •  .■* .  •  .VA  .  • . 
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In  Case  1,  Equation  E12  reduces  to 


/3  =  0 


and  the  roots  are 


Yl  =  Y2  =  Y3  *  0 


Therefore  Equations  E45  yield 


x2  _  x2  _  .2  _  _  a 

I  "  *11  '  AIII  "  3 


so  that 


X21  =  X31 


and  Equations  E45  yield 


Se 


i$ 


where 


Case  2a:  r  -  3  ?  0 
For  Case  2a,  we  set 

X11  =  /r  '  3 

X12  =  ~Jr  ‘  3 

r 

X  =  /§ei2  =  ^  cos!  +  i  sin- 


-KSe 


22 


=  -KSjcos|  +  i  sin| 


$ 


Si  ■  ^e‘i2  ■  ^ 


$  •  .  $ 
cos|  -  1  sin| 


x32  =  Vse~i2  =  /s 


<6  .  *.  I 

-cos|  +  i  sin; 


and  then  set 


X2  =  ✓§  cos| 

X3  «  /S  sin| 

so  that  the  mode  functions  are 

=  cosh  X^9 
W2  =  sinh  X^Q 
W3  =  cosh  X20  cos  X30 
W4  =  cosh  \2q  sin  *3® 

W5  =  sinh  X20  cos  X30 
Wg  *  sinh  X20  sin  X30 


so  that  the  mode  functions  are 

W|  =  cos  Xj9 

[F43] 

W2  =  sin  X]0 

[F44] 

W3' «  cosh  \2e  cos 

[F36  bis] 

W4  =  cosh  X20  sin  X30 

t F 37  bis] 

W5  =  sinh  X20  cos  X30 

[F38  bis] 

Wg  =  sinh  X20  sin  X30 

•  * 

[F39  bis] 

Case  3:  P  =  0;  Q  <  0 

In  Case  3,  Equation  E12  reduces  to 

y3  +  |  Q  |  =  0 

[F45  ] 

so  that 

Y3  =  -  l  Q  1 

t  F  46  ] 

Now  let 

r  *  |  Q  |  1/3  >  0 

[F16  bis] 

Then 

yi  -  -r 

[F47] 

TT 

n3  rf  r- 

V2  -  re  -  -  1  +  i/3 

[F48] 

-i- 

y3  =  re  3  =  |[l  -  i/a] 

[F49] 

and  Equations  E45  yield 

where 


We  now 


W.TW.iRVmn 7T TOT "■ 7 ’■»*>' » V» ir» v :-t:v.wjiwrjii'TOr>YJT.vvl v v v 

</ 


VII 


r  al  ,  r/3  _  i  4> 

2  ~  3J  +  1_ 2~  =  Se 


III 


-  [i  - 1]  -  4  ■  -■'* 


s2  = 


r  a] 2  3r^  2  ra  a^ 

2  ”  3  4  ~  r  ”  3  9 


<(>  «  cos 


-1 


r  a 
2  '  3 


set 


11 


•’/ 


r  +  3 


X, _  •  -i/r  +  ; 
12  V  3 


$ 


X,,  =  /§ei2  =  /§  cosf  +  i  sinf 


21 


$ 


X22  =  "*^el2  = 


cos 


2  +  1  Sinf] 


X 

X  *  /§e  i2  =  /sfcosf  -  i  sinf] 


[F51] 


[F52] 


[F53] 


[F54] 


[F55] 


[F56] 


[F27  bis] 


[F28  bis] 


f! 

S 


V 

N 

0 

:5 


6 


/■ 


t 


£ 


[F29  bis] 


X32  =  '^e  a  ^S[-cosf  ♦  1  sinfj 


[F30  bis] 


and  then  set 


'  a 
r  +  3 


X2  =  /S  cos^ 


X3  .  Vs  sin| 


so  that  the  mode  functions  are 

=  cos  Xj0 
W2  =  sin  X^fl 
W3  «  cosh  X20  cos  X30 
W4  *  cosh  X20  sin  X30 
W5  =  sinh  X20  cos  X30 
Wg  =  sinh  X20  sin  X30 

Case  4:  P  >  0;  Q  =  0 

In  Case  4,  Equation  E12  reduces  to 


[F57] 


[F32  bis] 


[F33  bis] 


[F43  bis] 
[F44  bis] 
[F36  bis] 
[F37  bis] 
[F38  bis] 
[F39  bis] 


Now  let 


y3  -  My  =  y[y2  -  lpU  *  0 


-  m 


yi  =  0 
y2  «  r 

y3  =  -r 


[F  58  ] 


[F59] 


[F60] 

[F61] 

[F62] 


•  v  -.1-  •••' 


x3i  -  Vr  +  5 


[F70  bis] 


X32  3  -Vr  +  3 


[F71  bis] 


and  then  set 


[F72  bis] 


X2  =  /  3  ~  r 


tF81] 


X  *  /  r  ♦  2 
3  V  3 


[F74  bis] 


so  that  the  mode  functions  are 


*1  *  cos  x^e 
W2  =  sin  X|0 

«3  =  COS  X20 
W4  =  sin  X20 
=  COS  X30 
Wg  =  sin  X30 


[F43  bis] 
[F44  bis] 
[F82] 
t  F83  ] 
[F77  bis] 
[F78  bis] 


Case  5:  P  <  0;  Q  =  0 

In  Case  5,  Equation  E12  reduces  to 


y3  +  1 p I y  =  y[y2  +  I p I ]  =  0 


[F84] 


a22  =  -vi>e  =  -yb|cos2  +  i  sm2j 


X^  =  /Se  i2  =  /s  cos^  -  i  sin^j 


X32  =  -t^Se  12  =  l/s  ~cos2  +  i  sin2 


and  then  set 


Wl 


X2  =  COS2 


X3  .  ✓§  sin2 


so  that  the  mode  functions  are 

=  cos  Xjfl 
W2  =  sin  X}0 
*3  =  cosh  \2d  cos  X30 
W4  =  cosh  X2d  si0  X30 
W5  =  sinh  X2©  cos  X30 
Wg  =  sinh  A20  sin  X30 

Case  6:  P  <  0;  Q  *  0 


In  Case  6, 


[F28  bis] 


[F29  bis] 


[F30  bis] 


[F72  bis] 


[F32  bis] 


[F33  bis] 


[F43  bis] 
[F44  bis] 
[F36  bis] 
[F37  bis] 
[F38  bis] 
[F39  bis] 


[F91] 


so  that  Equation  E20  yields 


^/r7 if 


and  therefore  Equations  E21  and  E22  yield 


A3  =  9  +  /r  >  0 


B3  =  9  -  /r  <  0 


Now  let 


fQ 


12 


+  KR 


1/3 


0 

8  =  |  -  VR 


1/3 


and  then  set 


Aj  =  A 


■  2n 


A2  =  Ae  3  =  |(-1  +  i/3) 


■  2n 


A3  =  Ae  =  "  f (1  ♦  i/3) 


B!  =  "  8 


i? 


B2  =  Be  3  =  |(1  +  i/3) 


-  i- 

B3  =  Be  3  =  |(1  -  i/3) 


75 


[F92] 


[F93] 


[F94] 


[F95] 


[F96] 


[F97] 


[F98] 


[F99] 

[F100] 


[F101] 


.  %  \  v  _ 

V~- 


/»><•> 


»v 

yyv 


ft  -  - 

:# 

^5 

ft;  . 


v;..\ 


*>  . 


•vs 


V-V 


'V.\ 
-  *  -  A 


.'•  .•> 


W‘ 

,w 


%-V'v 

vv« 

V^y- 

8$ 

1  «>’■ . 


-.  •.  •> 
•/sV 


.■>y. 


[  F 1 0  2  ] 


Then  Equation  £32  yields 


yl  =  A1  +  B1  =  A  -  8 


y2  =  A2  +  B2  =  2 


( -A  +  B)  +  i/3(A+B) 


y3  =  A3  +  B3  =  2 


(-A  ♦  B)  -  i/3(A  +  B) ] 


so  that  Equations  E45  yield 


»j  ■  (A  -  B)  -  | 


where 


II 


-A  +  B  _  a 
2  ".3 


./3  .  .  _  id) 

♦  i  ^ ( A  +  B )  =  Se 


'III 


-A  +  B  a 
2  '3 


i^f<A  +  B)  =  Se-1^ 


„2  A2  -2AB  +  B2  (-A  +  B)a  a2  3. .2 

S  = - : - -  - - —  +  —  +  r(A 


9  4 


(A 


♦  ab  ♦  b2)  ♦  ♦  £ 


<f>  =  cos 


-1 


A  -  B  a 
2  *  3 


Case  6a:  (A  -  8)  -  3  »  0 


+  2AB  +  B2) 


For  Case  5a,  we  set 


V 


X21  =  /Sel2  =  /s  cos|  ♦  i  Sin| 


[F27  bis] 


22 


$ 

=  -/§e’2  =  -/§  cos|  +  i  sin|j 


31 


4  r  1 

-  ✓se"’2  =  /s|cos|  -  i  sin^j 


•  -*■" *  ■  * 


-cos|  +  i  sin|j 


and  then  set 


A2  =  VS  cos| 
a3  .  K  sinf 

so  that  the  mode  functions  are 

=  cosh 
W2  =  sinh 

W3  *  cosh  \20  cos  X30 
W4  «  cosh  X20  sin  X30 
W5  =  sinh  X20  cos  X30 
Wg  =  sinh  X2®  sin  X30 


[F28  bis] 


[F29  bis] 


[F30  bis] 


[F113] 


[F32  bis] 


[F33  bis] 


[F34  bis] 
[F35  bis] 
[F36  bis] 
[F37  bis] 
[F38  bis] 
[F39  bis] 


[F114] 


Case  6b:  (A  -  B)  -  3  <  0 


For  Case  6b,  we  set 


X11  ■  Vf  -  <A  -  B> 


*12  =  'Vf  -  (A  -  B) 


*21  =  /i  cos  |  +  i  sin| 


A22  =  Vs  c os|  +  i  sin| 
A31  =  /S  cos|  -  i  sin^ 
X32  =  »/s  -cos|  +  i  sin2 


and  then  set 


3  -  (A  -  B) 


A2  =  ✓§  cos| 


A3  =  /S  sin^ 


so  that  the  mode  functions  are 

W|  =  cos  Ax0 
W2  =  sin  X^e 
W3  =  cosh  A26  cos  A36 
W4  =  cosh  A2©  sin  A30 


[F115] 


t F 2 7  bis] 


[F28  bis] 


[F29  bis] 


[F30  bis] 


[F116] 


[F32  bis] 


[F33  bis] 


[F43  bis] 
[F44  bis] 
[F36  bis] 
C F 37  bis] 


raw. 


KS5 


W5  =  sinh  X2$  cos  KjQ 
W6  =  si  nh  \20  sin  X30 


CF38  bis] 
[F39  bis] 


3 

[Ql 

I  2/3 

Case  7:  0  <  P  <  0 

2_ 

_ ;  Q  >  0 

In  Case  7, 


so  that  Equation  E20  yields 


[FI 17] 


[  F 1 1 8  ] 


and  therefore  Equations  E21  and  E22  yield 
A3  =  ^  +  /r  >  0 

B3  =  j  •  /R  >  0 


Now  let 


A  =  3  +  y r 


1/3 


8  -  [2  -  H 


1/3 


and  then  set 


A1  =  A 


i  — 

A2  =  Ae  3  =  ^(-1  +  i/3) 


[  F 1 1 9  ] 

[F120] 

[F121] 

[F122 ] 

[F97  bis] 

[F98  bis] 
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A3  *  Ae 


'3  A 


-  “(I  +  iV3) 


m  r»  v*  w  l>  /vvv^'v/v}.^ 

m 


[F99  bis] 


•  >.v.\ 

;.-v\ 


B1  *  B 


b2  *  Se  3  =  |(-1  +  i/3) 


®3  =  Be  3  =  ‘  |(1  + 


[F123] 


[F124] 


[F125] 


2?& 

m 


Then  Equation  E32  yields 


Vl  =  Aj  ♦  Bj  =  A  ♦  B 


A  ♦  B ,  _  .  rz. 

y2  =  A2  +  b2  =  2  (-!  +  i^3) 


y3  *  *3  *  b3  -  -  H"6'1  *  i/5> 


so  that  Equations  E45  yield 


XI  *  'A  *  8)  -  | 


■  -[v-5  * !]  *  4*  *  •)  -  «•* 

xni  ■  -[¥*  ‘  3]  -  4*  *  •>  - 


where 


[F126] 


[ F 1 27  ] 


[  F 1 2  8  ] 


[F129] 


[F130] 


[F131] 


2  .  (A  ♦  B)  iA_t. 
b  "  4  3 


..  nl2  (A  + 
=  (A  +  B)  +  - — 


g2  3  2 

+  -  +  -(A  +  B)2 


[F132] 


v^;v' 


•W  M 
\  vjCu 

CvCwffl 

>.vVl 

>&■*/ 


ivxvl 


-  .« , 

.vLvM 


.V.-.'.V 


•  .%  > 

.’y.v.\ 


f.  --  ■  ■■r.  r.  s 


»  A 


so  that  the  mode  functions  are 


Wj  =  cosh  Aj0 

W2  =  sinh  A}0 

W3  =  cosh  X20  cos  X30 

W4  =  cosh  X20  sin  X30 

W5  =  sinh  X20  cos  X30 

Ws  ■  sinh  X20  sin  X30 

a 

Case  7b:  (A  +  8)  -  3  <  0 
For  Case  7b,  we  set 


'll 


</F 


(A  +  B) 


'12 


=  -i/f  -  (A  .  8) 


A21  *  t^s|cos  |  +  i  sin^J 


x22  =  -/S 


4>  $ 

cos*  +  t  sin^ 


X3l  *  ^ 


d>  .  m 
cos|  -  1  sm| 


X32  ■  1/5 


d>  .  .  d> 

-cos|  +  i  sm| 


and  then  set 


X 


1 


(A  +  B) 


[F34  bis] 
[F35  bis] 
[F36  bis] 
[F37  bis] 
[F38  bis] 
[F39  bis] 


[F137] 

[F138] 

[F27  bis] 

[F28  bis] 

[F29  bis] 

[F30  bis] 

[F139] 


[F32  bis] 


Aj  3  rS  sin* 


C F33  bis] 


In  Case  8, 

0  < 


[! 


so  that  Equation  £20  yields 
0  i  i/r  = 


<  I 


[F117  bis] 


[ F 1 1 8  bis] 


so  that 

the  mode  functions  are 

.*v.. 

*1  * 
W2  = 

cos  A^fl 
sin  AjO 

[F43  bis] 

vx 

*  V*V  " 

%  *v  V 

•  v  v 
■VV*. 

[F44  bis] 

k  .< 

W3  = 

cosh  A28  cos  A38 

[F36  bis] 

W4  = 

cosh  A29  sin  A38 

[F37  bis] 

w5  = 

sinh  A20  cos  A30 

[F38  bis] 

we  - 

sinh  A20  sin  A38 

[F39  bis] 

Case  8: 

rQ]2/3 

0  <  P  <  3 1 2 1  0  <  0 

■A'-y 

■  V  •>  ' 


0\'-\ 


s\;.- 

.V  .  *  . 


and  therefore  Equations  E21  and  E22  yield 


A3  *  ^  +  »^R  <  0 


3  0*- 

8  =  |  -  KR  <  0 


[F140] 


C  F 1 4 1  ] 


Now  let 


sfc  = 


+  +  J(A  +  B)2 


{A  *  B)2.  <A  1_B ) a  +  |_ 


<p  =  cos 


-1 


A  +  B  a 
2  ~  3 


We  then  set 


11 


=  i/(A  +  B) 


+  ! 


X12  =  -i/(A  +  B) 


a 

+  — 
3 


$ 


\21  =  /Se12  =  /S 


A  .  ♦ 

cos|  +  1  sin* 


$ 


X22  =  -/Sei2  =  -t^S  cos|  ♦  i  sin| 


4> 


A3i  =  /Se'i2  =  /S 


A  .  A 

cos|  -  1  sm| 


4> 


X32  =  -|/Se  12  =  /s 


-cos 


i + 1  sinf]' 


and  then  set 


‘,-/J 


A.  »  y  (A  ♦  B)  +  | 


A.  =  /§  cos? 


[F156] 

[F157] 

[ F 1 58  ] 

[F159] 

[F27  bis] 

[F28  bis] 

[F29  bis] 

[F30  bis] 

[FI 60 ] 

[F32  bis] 


[F33  bis] 


A3  =  /S  sin* 


so  that  the  mode  functions  are 


=  cos 

W2  =  sin  A^0 
W3  =  cosh  A26  cos  A30 
W4  =  cosh  A2O  sin  A38 
W5  =  sinh  A28  cos  A38 
Wg  =  sinh  A2a  sin  A30 


Case  9: 


P  > 


Q  *  0 


In  Case  9, 


so  that  Equation  E20  yields 


We  therefore  define 

u  *  hos~l 

and  then  set 

Uj  =  0) 


[F43 

bis] 

[F44 

bis] 

[F36 

bis] 

[F37 

bis] 

[F38 

bis] 

[F39 

bis] 

[F161] 


[F162] 


[F163] 


[E39  bis] 

[E40  bis] 


U)  ♦ 


[E41  bis] 


u 


3 


3 


and,  following  Equation  £44,  set 


y 


1 


cos 


[F164] 


[F165] 


y3  = 


[F166] 


The  following  diagram  shows  that,  since 
0  <  u)  $  \ 


C  F 1 67  ] 


and  then  set 


X1  =  /yi  '  3 


[F178] 


X  =  /y  -  2 
2  /  y2  3 


[F179] 


X3  =  /  3  “  V3 


[F180] 


so  that  the  mode  functions  are 


W}  =  cosh  \^9 
W2  =  sinh 
W3  =  cosh  X2@ 
W4  =  sinh  X20 
W5  =  cos  X30 
W6  =  sin  X30 


[F34  bis] 
[F35  bis] 
[F75  bis ] 
[F76  bis] 
[F77  bis] 
[F78  bis] 


_  >  0;  A*  <  0 

Case  9b:  I  'II 


In  Case  9b,  we  set 


'1  3 


[F172  bis] 


V  yl  3 


[F173  bis] 


vry2 


[Fi8i] 


-  V!  - 


[F182] 


=  v;  -  Vo 


[F176  bis] 
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X32  *  ’  '/ 3  Y3 

[ F 1 7 7  bis] 

and  then  set 

X1  =  /yl  "  3 

[F178  bis] 

CM 

> 

1 

II 

CM 

«< 

[F183] 

X3  55  y 3  ”  y3 

[F180  bis] 

so  that  the  mode  functions  are 

Wj  »  cosh  XjG 

[F34  bis] 

W2  *  sinh 

[F35  bis] 

W3  =  cos  \29 

[F82  bis] 

W4  *  sin  Xo0 

[F83  bis] 

W5  =  COS  A30 

[F77  bis] 

Wg  =  sin  X30 

[F78  bis] 

«* 


Case  9C : 


<  0 


In  Case  9c,  we  set 


> 

II 

f 

< 

*-» 

[F-184] 

X12  =  _1/ 3  "  V1 

[F185] 

X  a  i  A  -  y  . 

[F181  bis] 

a21  V  3  y2 

X22  =  _1/ 3  "  y2 

[F182  bis] 

90 


.  l.- _ i 


[F176  bis] 


X31  V  3  “  y3 

X32  =  ~'/ 3  “  V3 
and  then  set 

xi  ■  A  •  yi 

X2  3  A  -  V 

X3  =  l/ij  "  V3 

so  that  the  mode  functions  are 

W1  =  cos  AJ0 
W2  *  sin  A^0 
W3  =  cos  A20 
W4  =  sin  X20 
W5  =  cos  X30 
W6  =  sin  X3e 


[F177  bis] 


[F186] 


[F183  bis] 


[F180  bis] 


[F43  bis] 
[F44  bis] 
[F82  bis] 
[F83  bis] 
[F77  bis] 
t F7 8  bis] 


Cases  1-9  (Summary) 

For  convenience,  the  above  cases  are  summarized  below. 


Case 

1 

2 

2a 


Criteria 
P  =  Q  =  0 

P  =  0;  Q  >  0 

r  -  a/3  *  0 


2b 


r  -  a/3  <  0 


Criteria 
P  >  0;  Q  =  0 


r  -  a/3  *  0 
r  -  a/3  <  0 

P  <  0;  Q  =  0 
P  <  0;  Q  *  0 

(A  -  B)  -  a/3  *  0 
(A  -  B)  -  a/3  <  0 

[012/3 

0  <  ?  <  3[2]  '  *3  >  0 


(A  +  B)  -  a/3  0 

(A  +  B)  -  a/3  <  0 


0  <  P  <  3 


ral 

12 


2/3 


Q  <  0 


P  >  3 


[Q12/3 

UJ  ' 


Q  *  0 


A*  >  0;  A^  >  0 

Aj  A  0;  A^  <  0 


<  0 
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APPENDIX  G 


SYMMETRY 


Any  function,  f(x),  can  be  expressed  as  the  sum  of  a  symmetric  (even) 
tion,  f £ ( x ) ,  and  an  antisymmetric  (odd)  function,  fQ(x), 


f  (x)  =>  f£(x)  *  fQ(x) 


where 


fE(x) 

f0(x) 

From  Equations  G2  and  G3 
function . 


=  fjxj  ±  f.(-x) 

2 

-  f(x)  -  f(-x) 

2 

come  the  definitions  of  an  even  and  an  odd 


fE(-x)  =  f E ( x ) 

f  (-x)  =  -f  (x) 

o  o 


From  the  definition  of  a  derivative,  and  using  the  notation 


we  obtain 


L  = 


1  im 


Ax-o 


fE(-X)  =  L 


fe(-x)  -  f E ( -x-Ax) 


Ax 


=  L 


f£(x)  -  f£(x+Ax) 


Ax 


=  -  f£(x) 


fQ(-x)  -  L 


fQ(-x)  -  fQ(-x-Ax) 


Ax 


-f  (x)  ♦  f  (x+Ax)l 

— - - y -  -  t.  i  v,  \ 


func 

[Gl] 

[G2] 

[G3] 

[G4] 

[G5] 

[G6] 

[G7] 

mo  i 


APPENDIX  H 


MODE  SHAPE  DETERMINATION 


Equations  B5,  C3,  and  C4  express  the  centroidal  translation  and  rotation 
of  a  plane  cross  section  of  a  circular  arch  vibrating  in  a  single  mode  in 
separable  form, 


ul' 

X 

(u>  = 

U2 

=  T  (U)  =  T 

Y 

ri); 

Z 

[HI] 


where, 

U1 

u2 

r 

T 

X, 


referring  to  Figure  A1 

=  centroidal  tangential  displacement,  in  the  direction  of  increasing  0 
=  centroidal  inward  radial  displacement 
=  clockwise  rotation 

=  arch  cross-section  radius  of  gyration,  defined  in  Equation  A62 

=  modal  amplitude,  a  function  of  time  only 

Y,  Z  =  mode  shapes,  functions  of  the  arch  angle,  0,  only 


Equation  HI  is  the  general  definition  of  a  structural  vibratory  mode,  i.e.,  a 
condition  in  which  the  relative  displacements  of  all  points  on  the  structure 
remain  constant,  and  the  absolute  displacements  are  proportional  to  a  single 
scalar  function  of  time.  The  configuration  defined  is  that  of  the  arch  cross- 
section  centroidal  axis,  not  its  neutral  axis,  because  the  centroidal  axis 
location  is  independent  of  cross-section  rotation. 


When  the  arch  undergoes  free  vibration.  Equation  C15  shows  that  the  modal 
amplitude,  T,  is  harmonic  with  angular  frequency  p. 

T  +  p2T  =  0  [C15  bis] 

Appendix  F  shows  that  the  three  mode  shapes,  X  (for  centroidal  tangential 
displacement),  Y  (for  centroidal  radial  displacement),  and  Z  (for  plane  cross- 
section  rotation)  can  each  be  expressed  as  a  linear  combination  of  six  mode 
functions , 


where 


UJ)  =3x1  column  matrix  of  mode  shapes 
(W)  =6x1  column  matrix  of  mode  functions 
C  =3x6  rectangular  matrix  of  mode  shape  coefficients 

Each  of  the  six  mode  functions,  W-j  (i=l,6),  satisfies  Equation  E2,  which  is 
a  homogeneous,  linear,  ordinary  differential  equation  with  constant  coef¬ 
ficients,  of  order  six.  Equation  E2  is  derived  from  Equations  C18,  which 
are  three  independent,  homogeneous,  linear,  ordinary  differential  equations 
with  constant  coefficients,  of  order  two.  Equations  C18  are  the  fundamental 
equations  describing  arch  free  vibration,  and  Equation  E2  was  derived  from 
them  for  the  sake  of  convenience,  to  uncouple  the  system  (cf.  Equations 
C21-C23).  However,  the  price  of  uncoupling  the  system  is  that  Equation  H2  is 
too  general.  Unless  restrictions  are  placed  on  the  mode  shape  coefficients, 
Equations  C18  will  not  be  satisfied.  To  determine  what  restrictions  must  be 
placed  on  the  mode  shape  coefficients,  we  substitute  Equation  H2  into  Equation 
C18  and  obtain 

[H3] 


[C17  bis] 

[B16  bis] 


The  six  mode  functions,  W-j,  are  a  closed  set,  in  the  sense  that  differen¬ 
tiation  of  any  one  mode  function  yields  a  linear  combination  of  the  other 
five.  Equations  H3  are  thus  a  set  of  three  homogeneous  linear  equations  in 


transcendental  frequency  equation  generated  by  Equation  H7  depends  on  the 
frequency  parameter,  a.  Thus  in  a  trial  and  error  solution  to  find  the 
natural  frequencies  and  mode  shapes  of  a  circular  arch,  the  form  of  the 
frequency  equation  varies  with  the  trial  frequency  parameter,  a. 


APPENDIX  I 


ORTHOGONALITY  RELATIONS 


The  basic  equations  for  the  mode  shapes  of  a  circular  arch  are 
(L  ♦  <j2  i){u|  =  {o | 

where 


L 


A02-B 

(A+8)0 

-ED2+BF 


- (A+B)0  -E02+BF 
B02-A  -(BF+E)D 
(BF+E  )D  CD2-BF2 


[CIS  bis] 


[08  bis] 


and 


0 

A 

B 

C 

E 

F 


d_ 

de 


I 

k 

2{Uv) 


2/3C2 

R 


[A60  bis] 

[D3  bis] 

[D4  bis] 

[D5  bis] 
[D6  bis] 
[D7  bis] 


The  constants  Cj,  C2(  and  C3  are  functions  of  the  ratio  of  arch  depth  to 
radius,  as  shown  in  the  table  following  Equation  A58. 

To  investigate  the  modal  orthogonality  relations  we  evaluate  Equation  C18 
for  modes  i  and  j . 


(L  +  o.2I) 

(u|f 

-  fo} 

[11] 

(L  *  0  2I) 

-  10} 

ri2] 

Premultiplying  Equation  II  by  and  premultiplying  Equation  12  by 

{U}^  yields 

{U}^  L  {U} i  ♦  a.2{U}^{Uf  .  =  0  [13] 

<<  L  {\i\.  +  aj2{u|^{U}j  =  0  [14] 

Subtracting  Equation  14  from  Equation  13  yields 

{Uf]  L  {U|.  -  {u|J  L  (U(  .  ♦  [a.2  -  a.2  {u}{{u}  .  =  0 

or 

a2  -  a2j  {U|{{U}.  =  {U^  L  {uj  J  -  {u}^  L  {uj.  [15] 

It  remains  to  be  shown  that  the  RHS  of  Equation  15  is  the  derivative  of  a 
function  which  has  equal  values  at  0  =  0  and  0  =  Q,  so  that 


arid 


(U)T  L  <U>  .  -  (U)T  L  <U>  .  =  b’.  . 
i  ~  J  J  “  1  ij 


[16] 

[17] 


If  Equations  16  and  17  are  satisfied,  then  integration  of  Equation  15  on 
0  between  0  and  Q  yields 


0 


[ia] 


2  2 

When  i  *  j,  so  that  o.  #  , 


Equation  18  yields 


Q 

J  {U|J{U|.  d0  =  0  (i#j) 

o 


[19] 


Thus  the  i*h  and  j**1  modes  are  orthogonal  over  the  closed  interval 
0  i  Q  i  Q,  provided  Equations  16  and  17  are  satisfied. 


To  prove  Equation  16  we  expand  its  LHS  using  Equation  08,  and  obtain 

b'.  .  =  A  fx.X."  -  X  .X .  "]  -  (A  +  B)  f  X  .  Y  -  X.y!]  -  E  [x  .  Z  . "  -  X  .  Z  .  " 

ij  i  J  J  i  l  i  J  J  iJ  1  J  J  i 

♦  (A  +  B)  Y.x'.  -  Y.x'.  +  B  Y.Y.  -  Y.Y.  -  (BF  +  E)  Y .  Z 

I  i  J  J  ]  L  1  J  J  1  1  J 

-  E  * [z  .X  ."  -  Z.X."1  +  (BF  +.  E)  Iz.y'.  -  Z.y!]  +  C  [z.Z.”  -  Z  .Z 
.■'j  J1]  i  j  jiJ  L  i  j  j 

ii  ti  it  ii  1  ii  it 

=  A  X.X.  -  X.X.  +  8  Y.Y.  -  Y.Y.  +  C  Z.Z.  -  Z.Z. 

L  i  j  jiJ  L  1  j  J1!  i  j  jiJ 


-  (A  +  B)  X.Y.  +  X.Y.  -  X.Y.  +  X.Y. 

L  1  J  1  J]  1  J  1  J  i 


(BF  +  E)  Y.Z.  +  Y.Z.  -  Y.Z.  +  Y.Z. 

L  i  J  ’  Jj  l  J  i  J  1 

II  II  '  II  II  ' 

E  Z.  X.  -  Z.X.  -  Z.  X.  -  Z.X. 

1  1  J  1  J  J  l  J  1  J  ij 


X.X.  -  X.X.  =  X.X.  +  X.X.  -  X.  X.  -  X.X.  =  X.X.  -  X.X. 
1  Jj  1  1  J J  1  J  1  J  1  J  1  J  1  J  J  1 


Z.X.  -  Z.X.  =  Z.  X.  +  Z.X.  -  Z.X.  -  Z.X.  =  Z.  X.  -  Z.X. 
1  jj  l  1  Jj  1  J  1  J  1  J  1  J  1  J  ■  1  J 


Therefore,  Equation  110  can  be  written  in  the  form 


B.  .  =  A  X.X.  -  X.X.  +  B  Y.Y.  -  Y.Y. 
iJ  1  JJ  l  i  JJ  ’  Jj  l  ■>  J 


+  C  Z.Z.  -  Z.Z  .  -  (A  +  B)  X  .Y  .  -  X  .Y  . 

1  Jj  l  1  JJ  l  1  J J  l  J  i 


.  - .  -\ 


_•  *  *  •  «*>*••*■  .  *  * 


-  »  *  |[Vjj  -  [YjZjj  J 

♦E  [[vj]''-  [zixj )'  -  [zWi  *  hxil' 


CH3] 


Each  term  on  the  RHS  of  Equation  113  is  a  derivative,  so  Equation  16  is 
proven. 

.  • 

To  prove  Equation  17  we  expand  Equations  B8  and  B9. 


{F}  =  eVs  {5} 


{6}  =  i  (DI  -  GT)  {u> 


[B8  bis] 


[B9  bis] 


where,  to  avoid  confusion  with  symbols  already  used  in  this  appendix 


E*  =  Young's  elastic  modulus 


A*  =  arch  cross-sectional  area 


Substituting  Equations  B9  and  C3  into  Equation  B8  yields 


*  *  *  * 

(F|  =  -R*-  S  (DI  -  GT)  {u}  =  ^  (DS  -  SGT)  (u| 


*  If 

=  ~~  T  (DS  -  SGT)  !u} 


[114] 


-2/3C, 


k 

0/1 


[B2,B22  bis] 


sgt  = 


2{l+v) 


k  R 

2r(l+i>) 


[B13  bis] 


so  that 


T  k 

OS  -  SG  =  ttt - - 

—  —  2{l+u) 

-2V2C2D 


2(UV) 


-  2V^C2D 


k  R 

2r( 1+u) 


[115] 


Now  i  f 


=  (N}  T  .  N2  T 
No 


[116] 


then  substituting  Equations  115,  116,  and  C3  into  Equation  114  yields 


A(X  -  Y)  -  EZ 


B(X  +  Y  -  FZ) 


★  ★ 
E  A 

R 


E(Y  -  X  )  +  CZ 


By  rearranging  terms.  Equation  113  can  be  written  in  the  form 


[117] 


[118] 


[119] 


eu  -  mtj  -  yjI  -  ezjI  -  hi4h  -M  -  Ez< 


*-VVa/.VV4V'/  */V,V  ‘a'/.V,’/,'.*  .*  /,  - 


-\J 

■  vvvj 


i  .4 


$3 


v.'.' \« 

,r\ 

V  r 
•'V\* 1 

r. 

r>.V, 

/  V  V 


y.% 

„V.>/ 
•  V  V 
.v.V 

V  ^  , 

Vvs 


103 


[120] 


+  Y.B  X.  +  y'.  -  FZ.  -  Y.B  X.  +  Y*.  -  FZ. 

1  J  J  JJ]  l  J  l  1  1  ^ 

I  I 

+  Z.  E  Y.  -  x'J  +  Cz'.  -  Z.  E  Y.  -  X*.  +  CZ*. 

U  L  J  jJ  JJ  I  JL  l  1  .■>]  TJ 

Substituting  Equations  117,  118,  and  119  into  Equation  120  yields 

b’..  =  [{u|t{n|  .  -  {u}t{n}.]' 

1J  *  *  1 1*  ' J  ’  J  ’ij 


[121] 


and  substituting  Equation  121  into  Equation  18  yields 

Q  Q 

o*  "  CTj2]  J  <ull(u|j  de  =  {uiTfNjj  -  {ujT{M|.  =0  [122] 

o  E  A  ° 

The  conditions  under  which  Equation  122  is  valid  for  i  *  j  must  still  be 
determined . 

Consider  the  elastic  restraint  conditions 


{N}  =  Fkj  {U| 
{N}  =  -  Fxj  {U} 


(8  =  0) 

(6  =  Q) 


Substituting  Equations  123  into  Equation  122  yields 

a 

a.2  -  a.2j  J  {Uj^Ujj  d8  =  [{U}J  fKjfU^  -  {uj^  [«j  {U|  .]  =  0  [124] 


be 

,v.y 

i* 

V.’A 

$$ 
.  m 

[I23a] 

T'r.v 

.  VfV' 

[  1 23b  ] 

r 

Thus  Equations  123  are  the  conditions  under  which  Equation  19  is  valid. 

One  last  task  remains  to  complete  the  orthogonality  analysis.  That  is  to  find 
the  integral  of  the  inner  product  of  a  mode  shape  with  itself  over  the  closed 
interval  0  <  8  <  0.  When  i  =  j,  Equation  122  reduces  to  the  indeterminate 
form 


*  '  «  *,*  *„•  ~ 


0 


J 

f  {u|T{u!  d0  =  {ujT{Nf  -  {u|t{n| 

1  2E  A  <7 


[133] 


f  =  f  (x1e,  x2e,  x3e ) 


[134] 


«X-  ^ 

If  - 


J  |U|>|, 


d0  =  0 


[125] 


•  .v  .  *  . 


As  suggested  by  Rayleigh  (Refs.  13  and  14),  L’Hospital’s  Rule  can  be  used  to 
evaluate  the  integral  in  question  by  considering  a  set  of  functions  identical 
in  form  to  { U }  and  {n|,  except  with  a  the  independent  variable  instead  of  0. 
Thus  if  we  set 


v>’> 

VV/ 

y'/j 


w> 


a  .  =  a 
J 


a .  =  a  +  do 


V-«-> 

|U|j  =  |ul 


{U}.  =  { U (  +  {U|  da 
{N| j  =  {N} 

{N}.  *  {N|  +  {N|  da 

then  applying  L'Hospital's  Rule  to  Equation  122  yields 
Q 

2a  da  [  {U|T{U(  d0  =  ^  [{u}T{n{  -  [u|T{n|] 

'  F  4  -In 


[126] 

v>) 

[127] 

»_ _ 

:-*• 

[128] 

[129] 

i 

[130] 

_v 

y.v> 

[131] 

-S-  > 

17* 

[132] 

*•  .V 

Vv> 

vv. 

*  *4  <»  I 


.s.V 


The  derivatives  with  respect  to  a  on  the  RHS  of  Equation  133  can  be  evaluated 
as  follows.  All  elements  of  { U }  and  { N. J  have  the  general  form 


tr'"i 


The  orthogonality  condition,  142,  is  used  in  Appendix  J  to  obtain  transient 
vibration  solutions  by  expressing  both  the  arch  displacements  and  the  forcing 
functions  in  modal  series  form.  See  Equations  J2  and  J3. 


APPENDIX  J 


TRANSIENT  ANALYSIS 


Equations  B17  and  C7  describe  the  transient  response  of  a  circular  arch. 
The  basic  equation  is 


<u)  -  [f]  L  <u>  =  -L-ib) 
.|RJ  “  pA 


Both  the  displacements  and  the  externally  applied  loads  are  assumed  to  be 
expressible  in  a  modal  series  expansion  of  the  form 


<u>  =  r  T.(U) ■ 

L —  J  J 

j=l 

00 

(b>  =  y~  h  .  (u> . 

L —  J  J 

j=l 

where  the  mode  shapes  satisfy  the  relations 


HU).  =  -<7.{U)  .  =  -  r  i2  { U  >  .  (no  sum) 
~  J  J  J  £  j 


[C16] 


[  (U)T(U)  ,d9  =  6.  . 
J  i  J  ij 


[142] 


Premultiplying  Equation  J3  by  (UK  and  integrating  on  0  from  0  to  Q  yields 


H.  =  |  { U > T ( b ) < 


Substituting  Equations  J2,  J3,  and  C16  into  Equation  J1  yields 


00  00  00 

Y~  T.(U>.+  y  P2  T  ,{U)  .  =  ^7  Y~  H  .  (U)  . 

L—  J  j  pA  L —  j  j 

j=l  j=l  j=l 


Premultiplying  Equation  J5  by  (U)^  and  integrating  on  0  from  0  to  £ 


T.  ♦  p2T.  =  i-H. 
l  it  pA  i 


so  that,  assuming 


T.(0)  =  T.(0) 


0 


the  solution  is 


t 

Ti(t)  =  pA  j  Hi(T)  Sin  Pi(t  "  T)dT 
0 


[J5] 

yields 

[J6] 

[J7] 

tJ8] 


which  completes  the  formal  analysis. 


APPENDIX  K 


VISCOELASTIC  ANALYSIS 


es 


The  dynamic  analysis  of  a  viscoelastic  circular  arch  turns  out  to  be  a 
fairly  simple  extension  of  the  results  for  an  undamped  elastic  arch,  when 
the  rate-sensitive  behavior  is  governed  by  a  single  parameter. 

When  the  arch  is  viscoelastic,  Equations  A40  and  A41  have  the  forms 


*«  [  a  «  ] 

*n  •  Eei  *  °  ei  ■  £i£i  *  r'u 


**12  =  Gy12  *  ny12  =  G[y12  +  G  y12j 


If  we  assume  that 


then  Equations  K1  and  K2  can  be  written  in  the  form 


all  =  E  (ei  +  5ei> 


a12  =  G  (y12  +  ^Y12] 


and  therefore  Equation  B8  has  the  form 


{ F }  =  EAS  (6)  +  £<S) 


and  Equations  817  and  C7  have  the  form 


Equation  C6,  which  describes  free  vibration  in  a  normal  mode,  has  the  form 


T<U>  -  §  (T  +  ?t)L<U>  =  (0) 


and  Equations  C12  and  C14  have  the  forms 


Tim)  -  |J  (T  ♦  5t)^Uj-1L  (U)  =  (0) 


— ^ -  <m>  =  \  <U}  =  -p2<m) 

T  +  lRJ  J 


[K10] 


Equation  Cl 5  has  the  form 


T  +  +  p2T  =  0 


[Kll] 


but  Equation  C16  remains  the  same. 


L  (U)  * 


<U)  *  -<T<U) 


[C16  bis] 


Equations  114  and  116  have  the  forms 


(F)  =  S^~(T  +  tf>(°§  -  sgt)(u) 


[K12] 


<F>  .  (T  +  ?f)<N) 


[K13] 


Equations  123  remain  the  same,  and  substituting  Equations  123  into  Equation 
K13  yields 


<F)  *  (T  +  $t)^Kj(U>  =  j  (u)  +  *{u) 
(F)  =  -(T  +  ) ^K, (U>  =  - ^Ki  (u)  +  ?(u> 


(0  =  0)  [K14a] 


(0  =  Q)  [ K 1 4b ] 


Thus,  in  order  for  the  modal  orthogonality  conditions  to  remain  valid,  the 
support  restraints  must  also  be  viscoelastic,  with  £  as  the  single  rate  sen¬ 
sitivity  parameter.  This  will  guarantee  that  the  form  of  Equation  124  will 
be  preserved. 

The  transient  analysis  problem  is  approached  by  substituting  Equations  J2 
and  J3  into  Equation  K7,  which  yields 


Y~  f  .  { U }  .  +  Y-  p2T.  +  $t.<U).  =  V~  H  .  <U>  . 

t—  J  J  £ —  J  L  J  Jj  J  pA  4_  j  j 

j-1  >1  j=l 


[K15] 


Premultiplying  Equation  [K15]  by  { U >  and  integrating  on  0  from  0  to  Q  yields 


2°  2  1 
T.  +  £p^T .  +  Pn.=  ^tH. 
i  'll  pA  i 


[K16] 


The  form  of  the  convolution  integral  which  is  the  solution  to  Equation  K16 
depends  on  the  values  of  £  and  p-j.  When  ?Pi  <  2.0  the  modal  damping  is  less 
than  critical;  when  £p,  =  2.0  the  modal  damping  is  critical;  and  when 
^Pi  >  2.0  the  modal  damping  is  greater  than  critical. 


APPENDIX  L 


NUMERICAL  ANALYSIS 

The  most  difficult  computational  phase  of  the  arch  modal  analysis  is  finding 
the  roots  of  the  frequency  equation  (H7).  A  sample  problem  for  a  Timoshenko 
beam  was  solved  numerically  on  the  AFWL  CRAY.  The  solution  involved  eleven 
program  segments,  which  were  stored  in  six  files.  These  program  segments 
are  described  below.  The  function  used  was  Equation  5-14  of  Reference  2. 

Program  TEST 

This  is  the  main  program  in  the  solution.  It  performs  all  of  the  I/O, 
determines  the  interval  and  tolerance  for  root  finding,  and  specifies  the 
number  of  roots  sought . 

Function  PRPCOM 

This  function's  purpose  is  to  prepare  the  common  region  PARAM,  which  con¬ 
tains  the  function  parameters  C12,  C22,  R,  L,  and  THETA.  The  meanings  of 
these  variables  are  given  below,  along  with  the  other  variables  of  PRPCOM. 


C12 

c5 

square  of  longitudinal  wavespeed 

C22 

square  of  shear  wavespeed 

R 

r 

radius  of  gyration 

L 

L 

length  of  beam 

E 

E 

beam  modulus  of  elasticity 

G 

G 

beam  elastic  shear  modulus 

RHO 

P 

beam  density 

I 

I 

beam  moment  of  inertia 

KPRIME 

K' 

shear  deformation  coefficient 

R8IG 

R 

rotational  beam-end  restraint 

THETA 


9 


EI/RBIG 


Subroutine  ZEROIN 

This  is  a  modified  version  of  the  ZEROIN  routine  found  in  Reference  16.  The 
modifications  were,  necessary  to  allow  the  routine  to  search  for  more  than 
one  root.  The  roots  are  stored  in  the  array  ROOT  as  they  are  found,  and  the 
routine  continues  searching  until  it  lias  found  the  desired  number  of  roots, 
or  until  it  cannot  find  any  more.  For  a  discussion  of  the  methods  used  to 
find  the  roots,  see  Reference  16. 


I 


Function  SGN 

This  simple  function  returns  1.0  if  its  argument  is  positive,  0.0  if  the 
argument  is  zero,  and  -1.0  if  the  argument  is  negative.  This  is  used  by 
ZEROIN. 

Function  G 

This  routine  uses  deflation  to  allow  ZEROIN  to  find  more  than  one  root  of 
the  function.  Deflation  is  a  technique  by  which  previously  found  roots  are 
eliminated,  but  other  roots  remain.  The  basic  idea  is  to  find  the  roots  of 
g(x)  =  f  ( x ) /Tl  ( x  -  r.j),  where  r-j  indicates  a  previously  found  root,  and  TI 
denotes  a  product.  In  thisway,  the  roots  are  eliminated.  Note  that  extreme 
care  must  be  taken  when  evaluating  the  function  at  points  close  to  a  pre¬ 
vious  root. 

Function  FACT 

This  simple  function  returns  the  factorial  of  its  integer  argument.  It  is 
used  by  Function  G. 

Function  F 

This  is  the  function  whose  roots  are  of  interest.  The  function  parameters 
are  contained  in  the  common  region  PARAM.  For  the  trial  function  they  were 
used  to  find  XI  (£)  and  GAMMA  (y),  which  are  used  to  find  K1  and  K2 .  The 
function  is  then  evaluated  piece  by  piece,  since  it  is  quite  complicated. 

Function  FI 

This  function  is  used  to  find  XI  and  GAMMA. 


Function  F2 

This  function  is  used  to  find  K1  and  K2 . 

Ref.  16.  Forsythe,  G.  E.,  M.  A.  Malcolm,  and  C.  B. 

Mathematical  Computations,  Prentice-Hall 


oler,  Computer  Methods  for 
(1977).  ~ . . . 


Function  DERIV 

As  mentioned  above,  care  must  be  taken  when  evaluating  the  function  g ( x )  = 
f(x)/jr(x  -  r-j)  near  x  =  r-j,  a  previous  root.  At  such  a  point,  [.’Hospital's 
Rule  must  be  applied.  This  means  that  the  derivative  of  f(x)  at  x  -  r ,  must 
be  found.  This  function  calculates  the  derivative  of  f ( x )  at  some 
point . 

A  Closer  Look  at  G  and  DERIV 

« .0 

As  mentioned  above, 

g(x)  s  f (x)/7t(x  -  ri)  r_; 

where  f(x)  is  the  function  of  interest  and  r-j  is  a  previously  found  root. 
This  definition  allows  us  to  ignore  all  previous  roots  without  losing  any 
others.  A  simple  example  will  help  illustrate.  If 

f(x)  »  (x-1 ) (x-2) ( x— 3 )  [12] 

then  the  roots  are  x  =  1,  2,  3.  Suppose  that  in  our  first  search  we  find 
that  x  =  3  is  a  root.  Then 

g(x)  *  [ (x-1) (x-2) (x-3) ]/(x-3)  =  (x-1) (x-2)  [L3] 

We  have  eliminated  the  root  x  =  3,  but  have  retained  x  =  1  and  x  =  2,  so 
that  one  of  them  will  be  discovered.  Note  that  once  all  the  roots  have  been 
discovered,  g(x)  a  1,  so  no  more  roots  will  be  mistakenly  discovered. 

This  method  works  equally  well  for  functions  other  than  polynomials,  but  it 
is  more  difficult  to  demonstrate.  The  problem  with  the  method  is  that  when 
ZEROIN  tries  to  evaluate  the  function  at  a  previous  root,  it  obtains 

g ( x )  a  f  ( x ) / ( x— r  -j )  =  0/0  [L4] 

which  the  computer  cannot  evaluate.  However,  through  the  application  of 
[.'Hospital's  Rule,  we  obtain 


g(x)  a  f(x)/7t(x 


ri)  =  f '  (xj/f^oix-n)]]] 


*  r  ir  »■  .  »  •  --  ^  ~  ~r  -r  ^  _•»  -r  -n  _r  ^r-  ^  j-  .- 


The  method  used  to  compute  derivatives  is  quite  simple.  It  is  based  on  the 
slope  interpretation  of  the  derivative: 


~  I 

dx  |x=xQ 


f(xQ  +  Ax)  -  f ( Xg  -  Ax) 
2Ax 


[L6] 


or 

.  f(x1)  -  f(x1  -  h) 
a  _ 

where  x^  =  xQ  +  Ax,  h  *  2Ax  . 


[L7] 


Higher  order  derivatives  are  found  simply  from  lower  order  derivatives 


df(n-D 

d,(n-l) 

df'"1 

• 

[d<n-l>x 

X  a  Xx  d(n_1)X 

x  =  x^  -  h 

d(nl* 

X  =  xn 

h 

[L8] 


This  procedure  is  implemented  easily  using  the  array  DFNDNX.  In  the  first 
step  of  the  routine,  enough  points  (n  +  1)  are  chosen  to  evaluate  the  nth 
derivative.  The  points  are  evenly  spaced  (h  is  the  interval  width  for  adja¬ 
cent  points),  and  the  point  at  which  the  nth  derivative  is  desired  is  at  the 
midpoint  of  the  interval.  As  the  points  are  chosen,  the  function  is  eval¬ 
uated  at  each  point;  and  these  values  are  stored  in  DFNDNX (l,i),  i  =  1, 
n  ♦  1.  In  the  third  step,  each  adjacent  pair  of  values  is  used  to  compute  a 
single  value  for  the  derivative.  This  step  is  repeated  until  there  is  only 
one  value.  This  is  the  n**1  derivative. 


For  example,  if  f(x)  x^  ♦  x^  +  x,  and  the  4*^  derivative  of  f  is  desired  at 
x  *  1.0,  then  DFNONX  will  contain  (if  h  =  0.01) 


OFNONX 


1 

2 

3 

4 

5 

f  ( .98) 

f ' (0.985) 

f  "  (0.99) 

f ’ ’ 1 (0.995) 

f iv(1.00) 

1 

2.8816 

5.8807 

7.9400 

6.0000 

0 . 0000 

f(-99) 

f ' (0.995) 

f ’ • (1.00) 

f"  ’  (1.005) 

■  •  •  • 

2 

2.9404 

5.9601 

8.0000 

6.0000 

f (1.00) 

f '  (1.005) 

f’ '(1.01) 

3 

3 . 0000 

6.0401 

8.0600 

f(1.01) 

f  ’  (1.015) 

•  •  •  • 

4 

3.0604 

6.1207 

f (1 .02) 

5  I 

! 

■ 

3.1216 

These  values  are  correct,  since 
f'(x)  =  3x2  +  2x  ♦  1 
f  "(x)  =  6x  +  2 
f  "  ’ (x)  *  6 
fiv(x)  =  0 


Values  for  less  well-behaved  functions  will  have  greater  error. 


Solving  for  Roots  of  Other  Functions 
When  seeking  the  roots  of  other  functions,  the  user  should  use  the  function 
f  to  compute  the  function.  Any  necessary  parameters  can  be  calculated  or 
assigned  in  PRPCOM,  and  passed  through  the  common  region  PARAM.  The  func¬ 
tion  f  must  be  a  function  of  one  variable.  The  user  must  be  aware  of  the 
meanings  of  the  arguments  of  ZEROIN. 


a  .the  lower  bound  of  the  interval 
b  the  upper  bound  of  the  interval 
f  the  function  of  interest 

tol  the  tolerance  (maximum  uncertainty)  of  the  independent 
variable  (root) 


**■»*«*< 
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nroot 

the 

number 

of 

roots 

desired 

root 

the 

roots 

iroot 

the 

number 

of 

roots 

found 

Thus,  in  order  to  solve  another  function,  both  f  and  PRPCOM  must  be  rewrit¬ 
ten,  as  must  lines  9-14  of- TEST.  All  other  code  is  completely  compatible. 
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real  function  derivCf ,n,x,h) 
implicit  non* 
real  f,x»h 
integer  n 

common  /machine  eps 

real  eps 

external  fndaps 

n  cannot  ba  graatar  than  10  unless  array  dimensions  ara  increased 

real  xlCll>,dfndnxCU,ll) 
real  tol , xtemp, dfn, dnx 
integer  i*j,maxj 

if ( aps. eq. 0.0) call  fndaps 

tol=2.0Xeps*x 

xl(l)=x-Cfloat(n)/2.0)*h 

xtemp-xl(l) 

dfndnxCl,  l)«f(xtamp) 

do  10  i=2»n+l 

xl(i)sxl(l)+floatCi-l)*h 

xtemp=xl (i ) 

df ndnxC 1 , i ) =f (xtemp) 

10  continue 

do  15  i=l,n+l 
15  continue 

do  30  i*2,n+l 
maxj  =n+2“i 
do  20  j^l.maxj 

df ndnx( i , j ) =df ndnxC i-1 , j+l )/h-df ndnxC i-1 ,  j  )/h 
20  continue 
30  continue 

do  50  i=l,n+l 
50  continue 

deriv=df ndnxC  n+1 » 1 ) 

return 

and 

subroutine  fndaps 

jf 

common  /machi n/  eps 
real  eps 

real  tol 

eps=l . 0 

10  eps=eps/2.0 
tol=l . 0+eps 
ifCtol.gt.l .0)goto  10 

return 

end 
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real  function  ftp) 
implicit  non* 
real  p 

common  /param/  cl2, c22, r, 1 , theta 
real  cl2,c22,r,l, theta 

real  fl,f2 
external  fl,f2 

real  xi , gamma , kl , k2 
real  fsign, temp 

real  xil, gl , col , co2, co3,co4,co5,co6 , tl , t2, t3, t4 
f=0.0 

if  (p.aq.O.O)goto  99 
fsign=-l . 0 
xi=flCfsign»p) 
fsign=+l . 0 
gamma-f 1 ( fsign,  p) 

fsign=-1.0 

terap=xi 

kl=f 2<  fsign, temp, p) 

fsign=+l . 0 

tempsgamma 

k2=f 2(  fsign, temp, p) 

xil=xixi 

if  ( xi 1. gt. 5000. 0)xi 1=5000.0 
gl=gammaXl 

col=2. 0 

co2=2. 0XthetaX(gamma+xiXk2/’kl ) 
co3=2 . 0*thetaX(xi+gammaXkl>'k2) 
co4=2 . OXthetaXthetaXgammaXxi 
co5=(k2/kl )X( l . 0+( ( thetaXxi )XX2) ) 
co6=(kl/k2)X( ( (thetaXgamma)xx2)-l .0) 

tl=colX(l-cosh(xil)Xcos(gl)) 
t2=co2Xcosh(xil)Xsin(gl) 
t3=co3x*inh(xil )»cosCgl) 
t4=(co4+co5+co6 )X*inh(xil )Xsin(gl) 

f=tl+t2-t3+t4 

continue 

return 

end 


real  function  fl(s,p) 


r 


real  s.p 

common  /par am/  cl2. c22. r, 1 , thata 
raal  cl2,c22,r,l, thata 

raal  tl,t2,t3,t4,t5 

tl=sX(l . 0/C12+1 . 0/c22)*p*p/2 . 0 

t2=(l . 0/c22-l . 0/cl2)*<p**4)/4. 0 

t3=C<p/r)XX2)/cl2 

t4=t2+t3 

t5=tl+sqrt(t4) 

fl=sqrt(t5) 

return 

end 


raal  function  f2(s,t,p) 
raal  s,t,P 

common  /param/  cl2, c22, r, 1 , thata 
raal  cl2, c22, r, 1 , thata 

raal  tl»t2,t3,t4 

tl*c22/<r*r) 

t2=tl-rho*rho 

t3=sxcl2*txt 

t4=tlxt 

f23<  t2+t3)/t4 

return 

and 
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raal  function  g(f,x.root,nroot) 

implicit  none 

raal  f,x 

raal  root(x) 

integer  nroot 

external  f 

external  fndeps 
common  /machin/  aps 

raal  aps 

raal  dariv 
integer  fact 
external  dariv, fact 


real  tol,h 
integer  i,n,nfact 
logical  atroot 

if  (eps.aq.0.0)  call  fndeps 
tol =abs( 2 . OXepsXx) 
if  (tol.lt.eps)  tol=eps 

atroot3. falsa. 
n=  0  • 

do  10  i3l,nroot 

if (abs(x-root(i) ) . gt . tol)  goto  10 
n=n+l 

atroot3 . trua. 

10  continue 

nfact=fact(n) 
h=0 . 001 

write(7, 100)nfact 
100  format(*  evaluating  limit* , i5) 

if  (atroot)  g=deriv(f,n,x,h)/float(nfactl 
if  ( .not. (atroot) )  g=f(x) 

do  40  i3l,nroot 

if  (abs(x-root(i ) ) .It. tol)  goto  40 
write(7,200) 

200  format(*  deflating  function*) 
g=g/(x-root(i) ) 

40  continue 

return 

end 

integer  function  faet(n) 
integer  n 

integer  i 

fact=l 
do  10  i=2»n 
fact=fact*i 
10  continue 

return 

end 
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raal  function  prpcom( dummy) 
implicit  none 
real  dummy 


II 


15  fa=g(f,a*root,iroot) 

fb=g(f . b, root , i root) 
e 

c  bagin  stop 
c 

20  c*a 

fc=fa 

d=b-a 

e=d 

30  if  (abs(fc) .ga.abs(fb))  goto  40 

a3b 
b=c 
csa 
fa=fb 
fb=fc 
fc=fa 


c 

c  convarganca  tast 

c  . 

40  toll=2. 0#eps*abs( b)+0 . 5*tol 

xm= . 5*(c-b) 

if  (abs(xm) . la. toll 1  goto  90 
if  (fb.eq.0.0)  goto  90 
e 

c  is  bisaetion  nacessary? 

if  (abs(e) . It . toll )  goto  70 
if  (abs(fa) . la.abs(fb) )  goto  70 
c 

c  is  quadratic  interpolation  possible? 
c 

if  (a.na.c)  goto  50 
c 

c  linear  interpolation 
c 

s=fb/fa 
p-2 . 0*xm*s 
q=l . Q-S 
goto  60 
c 

c  inverse  quadratic  interpolation 
c 

50  q=fa/fc 

r-fbyfc 
s=fb/fa 

p=s*( 2 . 0*xm*q*<q-r )-(b-a)x( r-1 . 0)  ) 
,  q*(q-1.0)»<r-1.0)*(s-1.0) 

c 

c  adjust  signs 
c 

60  if  (p.gt.0.0)q*-q 

p=abs(p) 
c 

c  is  interpolation  acceptable? 


if  <(2.0*p) .ge.(3.0*xm«q-abs(toll*q)))  goto  70 
if  ( p . ga . abs( 0 . 5Xe*q) )  goto  70 


ooo  o  o 


I 


d=p/q 
goto  80 

bisoction 

0  d=xm 
a=d 


complata  stop 


80  asb 

fa=fb 

if  (abs(d) . gt . toll )b=b+d 

if  (abs(d) . la. toll )b=b+sign( toll, x») 

fb=g( f ,b, root, i root) 

if  (CfbX(fc/abs(fc)l) .gt.0.0)  goto  20 
goto  30 
e 

c  dona 
c 

90  if  (sgnCfa)  .na.sgn(fb)  )goto  95 

if  (sgn(fc) .no.sgn(fb) )goto  95 
if.  (fb . eq . 0 . 0 )goto  95 
roturn 
c 

95  iroot=iroot+l 

rootCiroot)=b 
if  (iroot.aq.nroot)return 
c 

a=ax 
b-bx 
goto  15 
and 

raql  function  sgn(a) 
implicit  nona 
raal  a 

if  (a.eq.0.0)  goto  10 
sgn=a/abs(a) 
raturn 

10  sgn=0.0 

raturn 
and 
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program  tast 
implicit  nona 
raal  prp com 

axtarnal  f ,zaroin,prpcom 

raal  a,wprima, tol,root(100), dummy 

intagar  i , nroot, iroot 


II 


125 


opan(unit*7 , fila*' output* ) 
a  =  0 . 0 

wprim«=prpcom( dummy) 
tol=1.0a-6 

nroots20 

call  zaroin(a> wprima, f , tol » nr oot, root* i root) 
writa(7 , 90 )iroot 
do  10  i=l,iroot 

writa(7, 100)root(i) 

10  continua 

90  formate •  zaroin  found  '  *  i5*  '  roots') 

100  f ormatt al2 . 6 ) 
stop 
and 
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zaroin  found  8  roots 
0 . OQQOOOa+OO 
0.000000a+00 
0. 142336 a+05 
0 . 369554a+04 
0.147244a+04 
0 . 750657a+03 
0.202099a+03 
0 . 288620«+02 


